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Lecture  3 


GENERAL  DEFLAGRATIONS 

In  the  last  lecture  ve  examined  the  plane,  steady,  adiabatic,  premixed 
flame  and  deduced  an  explicit  formula  x'or  its  speed.  By  using  a  judicious 
choice  of  parameters  this  formula  can  be  made  to  agree  roughly  with 

f 

experiment;  precision  is  not  a  reasonable  goal,  given  the  crude  nature  of 
our  model.  Noteworthy  is  the  extreme  sensitivity  of  the  speed  to 
variations  in  the  flame  temperature:  an  0(l)  change  generates  an 
exponentially  large  change  in  flame  speed.  Such  variations  in  speed 
(caused,  for  example,  by  changes  in  mixture  strength)  are  not  excessive 
numerically  (at  least  for  fuels  burnt  in  air),  because  activation  energies 
and  fractional  changes  in  temperature  are  modest;  but  in  an  asymptotic 
analysis  they  present  a  potential  obstacle  to  discussion  of  multidimen¬ 
sional  and/or  unsteady  flames.  Then  significant  variations,  spatial  and/or 
temporal,  in  the  flame  temperature  can  be  expected  and,  if  the  sensitivity 
mentioned  above  is  any  guide,  there  will  be  correspondingly  large  spatial 
and/or  temporal  variations  in  the  flame  speed.  A  mathematical  framework 
ir.  which  to  accommodate  these  is  not  obvious.  (The  first  lecture  dealt 
with  special  circumstances  for  which  such  variations  were  manageable.) 

As  a  consequence,  attempts  to  discuss  general  deflagrations  have, 
for  the  most  part,  been  limited  to  situations  where  there  is  an  a  -priori 
guarantee  that  variations  in  the  flame  temperature  are  0(0  ^);  then  flame- 
soeed  changes  are  C(l)  and  present  r.o  mathematical  difficulties.  Two 
approaches  are  known  to  provide  the  guarantee  and  this  lecture  is  largely 
devoted  to  their  disclosure. 
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1.  The  Hydrodynamic  Limit 

At  the  end  of  section  2.4  the  steady  plane  deflagration  was  found  to 

Vave  a  thickness  5X/c  M  ,  and  this  may  be  taken  as  the  nominal  thickness 

?  r 

of  a  general  deflagration.  Ve  start  by  restricting  attention  to  waves 
whose  characteristic  length  (e.g.  minimum  radius  of  curvature)  is  large 
compared  to  their  nominal  thickness.  On  this  length  scale,  such  a  vave 

t 

is  simply  a  surface  across  which  jumps  in  temperature  and  density  occur 
subject  to  Charles's  lav  (as  in  appropriate  for  an  essentially  isobaric 
process) . 

If  the  ratio  of  the  two  scales  is  c,  then  on  either  side  of  the 
surface  the  appropriate  variables  are 

(x,y,z,t)  =  e(x,y,z,t);  (1) 

so  that  the  governing  equations  (2.13b, 19, 20)  becomes 

3p/3t  +  f-(cy)  =  0,  pDy/Dt  =  -?p  +  ePtV2*  +  |  7(?-y)],  (2) 

aDT/Dt  -  e?2T  =  0,  pDY/Dt-sI.-1  V  2  Y  =  0  (3) 


(We 

have  not  written  the  equations  for  components  other  than  the  single 
reactants  i  =  1,  and  the  subscript  1  has  been  dropped.)  As  e  •*  0,  ve 
have 


DT/Dt  =  DY/Dt  =0,  (4) 

i.e.  constant  values  of  T  and  Y  are  carried  by  the  fluid  particles. 

'.-.'e  conclude  that 


T 


(5) 


-3.3- 


everyv here  ahead  of  the  discontinuity  surface  if,  as  we  shall  suppose, 
these  constant  values  are  assumed  by  each  particle  at  its  point  of 
origin.  Likewise 

T  =  T.  ,  Y  =  0  (6) 

o 

everywhere  behind  the  discontinuity  since,  as  we  shall  see  presently,  these 

values  are  assumed  by  each  particle  as  it  leaves  the  flame.  Charles's 

law  (2.18a)  now  shows  that  p  has  the  constant  values  o ,»p.  =  p./a  on 

i  b  x 

the  two  sides  of  the  discontinuity,  where 


°  =  »f/p„  -  Vr  * 1  +  VTf 


(T) 


is  the  expansion  ratio  due  to  the  flame.  *ve  are  left  the  Buler's  equations 


2’V  =  0,  Pluf/?t  =  -7p 


(8) 


for  an  incompressible,  ideal  fluid,  i.e.  one  devoid  of  both  viscosity  and 
heat  conduction. 

The  two  ideal-fluid  regions  are  coupled  through  the  jump  conditions 


°f(vnf+V)  =  °b(vr.b+V)*  \f  =  yXb’ 

5f+pf(vnf+V)2  =  Vpb(vnb+V)2’  VYf  *  Tb; 


(9) 


(10) 


here  V  is  the  speed  of  the  deflagration  wave  back  along  its  normal 
(figure  1),  and  the  subscript  "X"  denotes  the  component  perpendicular 
to  ri,  i.e.  in  the  tangent  plane.  These  conditions  are  derived  in  the 
same  way  as  for  a  shock  wave  in  reactionless  gasdynamics,  i.e.  by 
integrating  the  basic  equations  (2.18b ,19,20)  through  the  flame.  Indeed, 
the  conditions  (9^  and  (10a)  are  identical  to  those  for  a  shock  since 
-hey  follow  from  the  same  continuity  and  momentum  equations.  The 
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requiremer.t  (lO'o)  can  also  be  recognized  as  a  Rankine-Hugoniot  condition, 
but  vith  kinetic  energy  neglected  and  a  heat-release  ter-,  (Y^.)  added. 

In  follows  from  the  combination 

pD(T+Y)/Dt  =  V2(T+L_1Y)  (11) 

of  the  basic  equations  (2.20). 

As  for  the  shock  vave,  these  jump  conditions  are  insufficient .  If 

the  state  f  immediately  ahead  of  the  vave  is  given,  there  are  five  scalar 

ecuatio.ns  for  the  six  unknown  scalars  p  (=p_T  / T  ),v  ,v  ,p  ,  and  V. 

b  i  f  b  r.b  jj3  b 

In  the  case  of  a  shock,  another  condition  is  imposed  from  outside  (such  as 
the  deflection  of  the  streamlines  at  a  sharp  body  or  the  pressure 
behind  the  vave  in  a  shock  tube).  Here  there  is  no  external  condition; 
ohe  deficiency  arises  from  discarding  information  by  using  only  the 
combination  (ll)  of  the  basic  equations  (2.20).  The  reaction  rate  then 
plays  no  role  in  the  derivation  of  the  jump  conditions.  Otherwise  stated, 
the  combustion  inside  the  vave  will  provide  information  about  the  burning 
rate  o _„V.r ,  i.e.  the  vave  speed 

*  “  \f  +  (12) 

Evaluation  of  VI  from  a  combustion  analysis  has  often  been  side¬ 
stepped,  even  though  it  should  be  considered  the  central  question  of 
premixed  flames.  Instead,  hypotheses  are  introduced;  the  simplest  is 
that  W  is  a  constant,  given  by  the  burning-rate  formula  (2.U3)  of 
steady,  plane  deflagrations.  This  hypothesis  is  justified  for  slowly 
varying  flames  (section  3)  when  L  =  1,  and  we  shall  use  it  in  lecture 
12.  3ut,  in  general,  it  is  not  acceptable  and  attempts  have  been  made 
(notably  by  Markstein)  to  modify  it,  in  particular  by  taking  into  account 
r.on-rlanar  characteristics  of  the  flame. 
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The  remainder  of  this  lecture  will  be  preparation  for  the  more  general 
combustion  analysis  that  follows  in  the  next. 

2.  Governing  Equations  for  the  Constant-Density  Approximation 

Although  the  formulation  can  be  carried  through  for  the  full  equations 
(2.18-20),  all  the  essential  features  are  preserved  under  the  assumption 
that  density  variations  due  to  the  presence  of  the  flame  are  negligible.' 

If  no  temperature  differences  are  imposed  on  the  flow,  the  velocity  field 
is  that  of  a  constant-density  fluid  and  can  be  calculated  in  advance;  ve 
shall  suppose  the  fluid  is  at  rest.  In  other  words,  we  shall  set 

p  =  1,  v  =  0  (13) 

in  the  full  equations  to  obtain 

3T/at  -  v2t  =  n-e-1^(T),  3Y/at  -  L_1v2y  =  -n  (l!0 

as  those  governing  the  combustion  field  under  the  constant-density  approxima¬ 
tion.  (All  equations  (2.20b)  except  the  first,  corresponding  to  the  single 
reactant,  can  be  omitted;  the  subscript  1  can  then  be  dropped.)  Note  that 
ve  have  added  a  term  -9  X'J'(T)  to  the  temperature  equation,  representing 
small  bulk  heat  loss. 

If  the  representative  mass  flux  M  is  chosen  to  be  the  burning  rate 
(2.1*3)  of  the  plane,  steady  (adiabatic)  deflagration,  then  the  reaction 
term  becomes 

a  /-n  2  2^  ^’o  L 

0  =  VYe  /x  with  V  =  Yf9  e  °/2LT^.  (15) 

Note  that  L  is  not  necessarily  equal  to  1  in  these  equations:  the  Lewis 
number  plays  a  very  important  role  in  the  analysis,  especially  for  unsteady 
flames.  Finally,  the  heat-loss  term  is  difficult  to  justify  in  a  multi- 
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diner.sior.al  context  (radiation  loss,  one  of  the  fev  legitimate  candidates, 
is  negligibly  small  unless  there  are  solid  particles  such  as  soot  in  the 
mixture);  but  for  quasi-plane  flames  it  can  represent  multidimensional 
effects  such  as  losses  to  sidewalls. 

We  shall  require  that 


T  -*  T  ,  Y  -*•  Yf  as  x  -*■  -® 


(16) 


ar.i  deal  exclusively  with  situations  where  equilibrium  prevails  behind  the 
flame  sheet,  i.e. 


Y  =  0  in  the  burnt  gas  (17) 

7r.e  temperature  behind  the  flame  will  be  close  to  the  adiabatic  flame 
temperature  (2.29). 

The  constant-density  approximation,  on  which  most  of  the  premixed  flame 
analysis  in  these  lectures  is  based,  clearly  provides  substantial  simplifi¬ 
cations.  It  can  be  Justified  as  a  formal  limit  in  which  the  heat  released 
by  the  reaction  becomes  vanishingly  small  (compared  to  the  existing  thermal 
energy  of  the  mixture).  Small  heat  release  can  be  due  to  either  a  scarcity 
of  reactant  (Yf  •*  0)  or  weak  combustion  (Tf  ■*  ®);  by  confining  ourselves 
to  dilute  mixtures,  we  have  already  assumed  the  former.  The  relevant 
parameter  is  the  expansion  ratio  (7);  asymptotic  expansions  in  o-l  provide 
a  formal  basis  for  the  approximation. 

Z.  Slow  Variations  with  Loss  of  Heat 

As  an  introduction  xo  the  more  complicated  analysis  of  multidimensional 
flames  consider  first  a  plane  flame  sheet  (figure  1),  looking  like  the 
adiabatic  deflagration  studied  in  section  2.U  but  moving  unsteadily  because 
cf  local  fluctuations  in  T  and  Y  (represented  by  the  time  derivatives). 
In  general,  the  flame  speed  can  be  defined  in  terms  of  the  mass  flux  of 


tie  mixture  through  the  sheet.  (This  is  a  veil-defined  concept  in  the  limit 
3  -*•  “i  for  9  finite  there  is  r.o  natural  definition,  except  vhen  the  combus¬ 
tion  field  is  steady  in  some  moving  frame.)  For  the  const ant -density 
approximation  adopted  here,  the  speed  is  Just 

V  =  -x*(t). 

"ote  that,  since  the  speed  is  not  defined  for  finite  9,  to  expand  it  in  the 
subsequent  asymptotic  analysis  vould  'oe  a  futile  gesture. 

Suppose  nov  that  changes  in  the  flame  speed  occur  on  an  0(9)  time 
scale,  i.e.  that 


t  =  t/3 


(19 


Ls  the  appropriate  (slow)  time  variable  to  describe  them.  Then,  for  an  ot 
r.ovir.g  vith  the  flame  sheet,  the  combustion  field  is  quasi-steady  to  leading 
:rder  (i.e.  steady  for  t  =  0(1)).  The  temporal  variations,  along  vith  heat 


Loss,  create  0(9  1)  perturbations,  and  hence  generate  only  0(9  1)  variations 


Ln  ohe  flame  temperature, 
'.•.her.  the  coordinate 


n  =  x-x*(t),  (20) 

based  on  the  flame  sheet,  is  introduced  the  basic  equations  (ll)  become 

9~‘3T/3t+Y  3T/3n-327/3r.2  +  9-1y  =  -9-13Y/3r  -V3Y/3n+t-132Y/3n2  =  fi.  (21) 

These  govern  the  motion  of  vhat  is  known  as  a  slowly  varying  flams  (SVF). 

To  integrate  the  equations  it  is  first  necessary  to  say  something  about  the 
flsme  temperature.  To  leading  order,  we  have 


V3(T+Y)/3r.  =  32(T+L-1Y)/3n2 


(22) 
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rrect  tc  leading  order.  To  the  sane  order  behind  the  flame  sheet,  T  is 
r.stant  (hence  shoving  that  the  derivative  in  the  result  (23)  is  zero,  as 
ticipated);  to  one  more  term  we  find 

T  =  7,  -  9 _1  [ V-~'o  ( 7.  ) n+T^6 ..  ( t )  ] ,  Y  =  0  for  n  -  0  (26) 

c  bo- 

writing  T  =  7,  in  the  5  ^  terms  of  the  temperature  equation.  Here  <}>#» 
oresentir.g  the  perturbed  flame  temperature,  is  as  yet  unknown. 

The  structure  problem  for  the  reaction  zone  determines  d*  as  a  function 
V.  This  problem  has  already  been  discussed  in  section  2.5,  where  the 
cression  (2.^7)  for  the  temperature  gradient  Just  ahead  of  the  flame 
eet  was  developed.  The  sar.e  gradient  can  be  calculated  from  ths  result 
ja),  leading  to  the  relation 

-d*/2 


-3.9- 


Clsarly,  there  is  the  same  temperature  sensitivity  as  for  steady  adiabatic 
deflagrations,  as  expected.  Moreover,  for  such  a  deflagration  the  perturba¬ 
tion  vanishes  and  7=1,  which  confirms  the  burning-rate  fornula 

v  -  •  O  )  • 

Another  relation  between  ^  and  V  comes  from  calculating  the 
change  in  enthalpy  of  the  mixture  between  its  fresh  and  burnt  states. 

For  that  purpose,  we  rewrite  equation  (22)  correct  to  0(e  h  before 


integrating  it  as  before,  to  obtain 

0+  -  . 
e_1  /  ~  (T+Y)dn  +  [ V(?+Y )  J  +  +  s”1 

J  T  — OO 


r  „  r  3T  .-1  3Y,0+  .  e. 

I  *  t  ^  *  t  ^1,  •  (28) 


The  integrals  can  be  evaluated  to  leading  order  by  means  of  the  formulas 
(25);  we  find 

/  ^-[Y^(eVn-eLVn)]dn  =  Y,V  /  n(eVn-eLVn)dn  =  -Y  (l-L_1)V-2V, 

•  o  T  x  I 


-v  x  l„-l 


/  y  (T_,+Y..e  n)an  =  [/  ii>(fT,+Y  e  )dv  ]V 

-  -  r\  - 


where  the  dot  is  used  to  signify  rate  of  change  on  the  t -scale.  The  formulas 
(2c)  and  boundary  conditions  (ic)  enable  the  remaing  terms  to  be  calculated; 


[V(7+Y)]^  =  V(T*-Tb)  =  -^'1TbV(()tf, 


r3T/an]0+  =  -0_1^(T  )V-1 ,  [3Y/3.n]0+  =  0. 

—  oo  D  —oo 

7::e  equation  (28),  in  which  all  terns  have  now  been  evaluated  to 


-nereiore  .3 1  ve s 


=  ?V_2-bV"3V  With  b  =  Y  (l-L-1)/T2; 


here 
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V  =  IC(T  )  +  /  i>(?r+?  e_V)dv]/?f\ 
0 


(30) 


the  tvo  terns  representin'  heat  lost  to  the  burnt  mixture  and  through 
the  side vails  ahead  of  the  flane  sheet,  -espectively . 

By  eliminating  d*  betveen  the  tvo  relations  (27)  and  (29) >  we 
obtain  an  equation  for  V,  namely 

bV  =  V3ZnV2  +YY.  (31) 


The  only  difference  when  the  constant-density  approximation  is  not  used  is 
a  more  complicated  formula  for  b.  The  crucial  property 

b  *  0  accordingly  as  L  *  1  (32) 

is  unaffected,  however.  Note  that  the  SVF  is  not  a  solution  of  the  general 
initial -value  problem  (only  the  value  of  V  may  be  prescribed  at  t  =  0); 
it  merely  describes  the  subsequent  behavior  of  any  flame  that  survives 
development  on  the  t-scale.  Thus,  a  prediction  of  instability  is  reliable 
but  r.ot  one  of  stability,  since  the  flame  may  have  already  lost  stability 
during  its  evolution  on  the  t-scale. 

Consider  first  the  steady  state  (figure  2)  determined  by  setting 

V  =  0  ir.  the  evolution  equation  (31),  i.e. 

V  =  0  or  V2  Zn  V2  +  Y  =  0.  (33) 

Cn  the  first  of  these  d#  =  +«,  so  that  the  perturbation  analysis  breaks 
down;  the  corresponding  nonunifornitv  has  never  been  treated.  The  second 
curve  provides  tvo  solution  branches  so  long  as  the  heat  loss  is  not  too 
l?.rge,  i.e.  H1  is  less  than  e  ^  (  =  0.363);  the  adiabatic  flame  speed 

V  =  1  is  attained  on  the  upper  branch  as  V  -*■  0,  so  that  it  is  plausible 
to  suppose  that  this  is  the  physically  relevant  one.  Ho  solution  exists 
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fer  ?  >  e  :  steady  combustion  car.r.ot  be  sustained  if  the  heat  loss  is 
too  large,  any  existing  flame  being  quenched.  It  is  interesting  that  the 
speed  of  the  flame  at  quenching,  namely  e  ‘  (  =  0.6C7)  times  its 
adiabatic  value,  is  completely  independent  of  the  nature  of  the  heat  loss, 
i.e.  the  form  of  the  function  y .  The  quenching  phenomenon  provides  a 
qualitative  explanation  of  the  Davy  safety  lamp:  the  vine  gauze  surrounding 
the  flame  is  an  effective  heat  sink,  preventing  the  propagation  of  the  flame 
beyond  its  confines. 

Equation  (31)  describes  the  evolution  of  plane  SIT s .  When  1  =  1, 

’o  is  zero  and  there  is  no  evolution:  equidiffusion  prevents  any  variation 
nr.  the  t -scale.  Ir.  fact,  since  the  equation  is  asymptotic,  there  is  no 
evolution  vhen  L  is  close  to  1,  i.e. 

r1  =  1  -S./e  with  l  =  0(1).  (3*0 

cut  then  a  treatment  on  the  t-seale  is  possible  in  certain  circumstances, 
leading  to  the  near-ecuidiffusion  flame  (DIF)  discussed  later. 

An  immediate  consequence  of  the  evolution  for  L  >  1  (b  >  0)  is  that 
the  flame  is  unstable:  any  deviation  of  V  from  its  value  on  the  upper 
branch  of  the  curve  ir.  figure  2  is  amplified.  The  same  conclusion  cannot 
be  dravn  for  L  <  1,  but  this  is  a  consequence  of  considering  planar 
disturbances  only.  Lecture  5  will  examine  the  linear  stability  of  plane 
deflagration  vaves  in  complete  detail,  and  find  that  plane  SVEs  are  unstable 
to  r.on-cianar  disturbances  for  L  >  1.  Thus,  the  SVFs  are  unstable  for  all 
values  of  L,  uhich  decreases  their  value  as  a  class  of  solutions  'but 
does  not  eliminate  them) . 

Mult i dimensional  Planes 

Consider  r.ov  situations  in  vhich  the  flame  sheet,  in  addition  to  being 
ur.soeady,  moves  ir.  a  non-planar  fashion.  The  goal  is  to  find  conditions 
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£0  tnan  vs  nay  writ 


=  fri  +  /  [vf(v+r-:)  -  If-ldn  -  e~-  /  ♦dn. 

** 1 0+  -=  ~L~  "  -CO 


.r.is  expresses  the  deviation  of  the  flame  temperature  T#  from  its  adiabatic 

value  in  terns  of  the  heat  lost  to  the  burnt  mixture,  the  transverse 

o 

diffusion  of  heat  and  reactant  up  to  the  flame  sheet ,  the  temporal  variations 
in  enthalpy  H  of  the  mixture  ahead  of  the  flame  sheet,  and  the  heat  loss 
up  to  the  flame  sheet . 

If  deviations  of  T#  from  are  to  be  0(0  ^),  the  right  side  of 

equation  (39)  aust  be  of  the  same  order.  This  is  guaranteed  when  the 
terms  in  3/3n,  V^,  and  3/3t  are  made  separately  small,  a  step  that  can 
be  taken  in  tvo  different  ways.  Cne  way  is  to  confine  attention  to 
disturbances  of  a  steady, plane  deflagration  that  vary  over  tines  and 
distances  0(0),  These  SV?s  are  a  generalization  of  the  ones  introduced 
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It  should  be  emphasized  that  S'.Ts  ar.i  ITEFs  are  restricted  classes  o 
solutions,  identified  by  setting  down  sufficient  (but  not  necessary) 
conditions  for  the  flame-temperature  variations  to  be  0(9  ■*■),  itself  a 
sufficient  condition  for  the  efficacy  of  our  asymptotic  method.  While 
these  classes  m  ay  be  the  only  general  ones,  special  circumstances  make  i 
possible  to  treat  other  prefixed  flames,  lack  of  time  prevents  our 
discussing  the  most  important  of  these,  namely  the  spherical  (premixed) 
flare:  symmetry  ensures  that  the  temperature  does  not  vary  at  all  over 
its  flame  sheet,  so  that  it  need  not  be  either  an  SVF  or  an  i!EF.  ('lever 
theless,  for  certain  parameter  values  it  is  an  SVF  and  for  others  an  HE” 
In  the  next  lecture,  the  equations  governing  the  SV7  and  the  NEF 
vill  be  derived  and  then  solved  for  a  basic  non-uniform  velocity  field: 
stagnation-point  flow. 
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Figure  Captions 


3.1  Notation  for  flame  as  hydrodynamic  discontinuity. 

3.2  Steady  flame  speed  V  versus  heat-loss  parameter  ¥  . 
Arrows  show  direction  in  which  speed  changes  for  L  >  1 
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